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Abstract 

We report the results of a search for the decay modes J/ip and B'^ J/ipD^ir^. The 

analysis is based on 140 fb~^ of data accumulated by the Belle detector at the KEKB asymmetric- 
energy e'*“e“ collider. No significant signals are observed and we determine the branching fraction 
upper limits B{B^ —> J/ipD^) < 2.0 x 10“® and B{B^ J/ipD^Ti^) < 2.5 x 10“^ at 90% 
confidence level. These results rule out the explanation of the excess in the low momentum region 
of the inclusive J/ip spectrum as intrinsic charm content in the B meson. The branching fractions 
of the corresponding nonresonant decay channels are also reported. 

PACS numbers: 13.25.Hw,14.40.Nd 


*on leave from Nova Gorica Polytechnic, Nova Gorica 
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The inclusive spectrum oi B ^ J/'ip X has been studied extensively and is consistent 
with the prediction of nonrelativistic QCD calculations except for an excess in the low 
momentum region Hi This momentum region corresponds to the J/'ip meson recoiling 
against particle systems with an invariant mass of ~ 2 GeV/c^. The observed excess below 
0.8 GeV/c corresponds to a branching fraction of a few times 10““^. 


Several hypotheses 
modes proposed in Ref. 



have been pr^osed to explain this excess. One of the deca’ 


L prOT 

J/V^Ap U, has been studied by BaBar [S| and Belle 


icay 

i 


The measured branching fraction, of order 10“®, is too small to account for the excess. 

Ghang and Hon proposed intrinsic charm (cc) in the B meson as an explanation. The 
intrinsic charm pair transforms into a cc hnal state when the B meson decays. The most 
promising decay modes are B^{dbcc) —> and B^{uhcc) —>■ J/%IjD^'k^. According 

to this hypothesis, if the intrinsic charm content of the B is not much less than 1%, the 
branching fractions of the above decay modes could be ~ 10“^, while the prediction from 
the standard QGD framework is of order 10“® @]. 

In this paper, we report on a search for the decay modes B^ —^ J/'ipD^ and B^ 
J/ipD^Ti^. The analysis is based on a data sample of 140 fb“^, which contains 152 xlO® 
BB pairs, accumulated at the T(4S') resonance with the Belle detector Q at the KERB 8 


GeV e and 3.5 GeV e"*" asymmetric collider [11 


ratio of the second to zeroth Fox-Wolfram moments 
The selection criteria for Jjjh mesons decaying to 
cal to those used in Ref. 



The Belle detector is a large-solid-angle magnetic spectrometer that consists of a three- 
layer silicon vertex detector (SVD), a 50-layer central drift chamber (GDG), an array of 
aerogel threshold Cerenkov counters (AGG), a barrel-like arrangement of time-of-flight scin¬ 
tillation counters (TOF), and an electromagnetic calorimeter comprised of GsI(Tl) crystals 
(EGL). These detectors are located inside a superconducting solenoid coil that provides a 
1.5 T magnetic field. An iron flux-return located outside of the coil is instrumented to detect 
Kl mesons and to identify muons (KLM). 

Events are required to pass the hadronic event selection criteria Q. To suppress contin¬ 
uum backgrounds (e+e“ —>• qq, where q = u, d, s, c), we require i ?2 < 0.5, where i ?2 is the 

m 

(where I = e, fi) are identi- 
To remove charged particle tracks that do not come 
from the interaction region, we require that the leptons originate from within 5 cm of the 
nominal interaction point (IP) along the beam direction. Both tracks are required to be 
positively identihed as leptons. In order to reduce the effect of bremsstrahlung or hnal 
state radiation, clusters detected in the EGL within 0.05 radians of the original e~ or 
direction are added in the invariant mass calculation. The J/'ip candidate is required to 
satisfy an asymmetric invariant mass requirement that takes account of the radiative tail: 
—150(—60) < Me+e-( 7 )(M^+^-) — mj/^ < 36(36) MeV/c^, where mji^ is the nominal J/'ip 
mass [l^. In order to improve the momentum resolution, vertex and mass constrained hts 
are then applied to the J/'ip candidates that pass the above selection criteria. 

For charged pion and kaon identihcation, the specihc ionization (dE/dx) in the GDG, 
the hight time measured in the TOF, and the response of the AGG are combined into a 
likelihood L^, where h stands for the hadron type (vr, K, p). A track is labeled as a kaon 
if Lk/{Lk + LtP) > 0.5 or a pion if L^/(L^ -|- Lk) > 0.3; the respective efficiencies are 90% 
and 92%, while the respective n/K misidentihcation rates are 10% and 13%. Tracks in the 
kaon sample with Lp/{Lp -|- Lk) > 0.99 are reclassihed as protons and thereby discarded. 
All tracks positively identihed as electrons are rejected. 

A meson candidate is reconstructed from a K~ and a tt’*' meson pair. The two 
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tracks must satisfy dr < 0.3 cm and \dz\ < 5 cm, where dr (dz) is the impact parameter 
perpendicular to (along) the beam direction with respect to the IP, determined run-by-run. 
We select candidates for further analysis within the mass window \Mk-k — '^z)o| < 150 
MeV/c^, where is the nominal mass Q. We also apply a vertex constrained £t to 
the candidates. 

For the charged pion in —>• J^ we apply the requirements with looser dr < 0.6 

cm, \dz\ < 5 cm, and tighter -|- L^) > 0.9 due to its low momentum. 

B mesons are reconstructed by combining a, J/xjj and a candidate for > J/ip 
and an additional pion with the same charge as the kaon from for B^ —>■ J/'ipD^TT^. 
To reduce combinatorial background, we impose a requirement on the quality (y^) of the 
vertex £t for the leptons from J/ip and the trajectory (and the tt’*' for the B^ case). The 
vertexing requirement retains 94% (78%) of the B^ {B'^) signal. 

We reject B^ —> 'ip{2S)K~^ [t/’(2*S') —> decay by requiring — 

Mi+i-{'y) to be outside of the ±15 MeV/c^ window around the nominal mass difference 
between 'ip{2S) and J/ip. We require | cos^^l < 0.8 to further suppress combinatorial back¬ 
ground, where 9b is the angle between the B flight direction and positron beam direction 
in the center-of-mass (cms) frame. 

We select B candidates by requiring that the beam-energy constrained mass (Mbc = 
F^b ) and the the mass difference (^kfVfb — ]\dB ^^^n) Q lie within the region 
Mbc > 5.2 GeV/c^ and —0.3 < AM^ < 0.2 GeV/c^; here, i^beam and Pb are the beam energy 
and B momentum in the cms, while Mb and rriB are the reconstructed and the nominal mass 
of the B meson. The signal region for B^ {B^) candidates is dehned as 5.27 < Mbc < 5.29 
GeV/c^, |AMs| < 19.0 (17.1) MeV/c^, \MK-K — mBo\ < 16.5 MeV/c^, corresponding to three 
standard deviation windows based on Monte Garlo (MG) simulation and the data control 
samples described later. The candidates outside of the signal region are used to determine 
the background components in the £t described below. 

After the selection, around 3.6% (42%) of B^ {B^) candidate events have more than one 
B candidate. In multiple candidate cases, we select the candidate with the best vertex £t 
quality. 

Background is divided into categories that have distinct shapes in the (Mbc, AM^, Mk-k) 
distributions: nonresonant background, two types of combinatorial backgrounds, and peak¬ 
ing backgrounds. In the case of nonresonant background, a B meson decays to the signal 
hnal state, but the and n~ mesons do not come from a decay. In the case of combi¬ 
natorial background, the reconstructed J/tp and [and a pion, for a B^ candidate] come 
from different B mesons (90% of the time) or from continuum events (10%). No peak ap¬ 
pears in the (Mbc, AM^) distribution; however, one subclass—cmb(DO)—has a peak in the 
Mk-k distribution (77° correctly reconstructed), while the other—cmb(fake DO)—does not 
(fake D^). 

The peaking background shows an enhancement in the Mbc signal region. For the B^ 
J/ipD^ signal, one source of this background is from B~^ —>■ 'ip{2S){xc)K~^ decay, where 
the J/pj and are combined with a pion from the B~ meson. Another source is from 
B^ —> J/'ipK^TT~Ti^ {B^ —>■ J/'ipK^Ti~TT^) decay where the second pion is missed. For the 
B^ —>• J/ipD^TT^ signal, one source of peaking background comes from the aforementioned 
B~^{B^) decay with the (tt^) replaced by a charged pion from the other B meson. Another 
is from the B^ —> J/'il)K~^7r~ decay combined with a from the accompanying B. For the 
B^ {B~^) signal, the hrst-mentioned peaking background distributes broadly around the 
AMb signal region, while the second exhibits a narrow peak shifted to negative (positive) 
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values by a pion mass. The narrow peak is excluded by limiting the fit region to AMb > 
-0.12 (< +0.12) GeV/cT 

The yields are extracted by maximizing the three-dimensional (3D) extended likelihood 
function, 


C 


e 

N\ 


N 


n 


^ PkiMU.AMh.M},,) 




where N is the total number of candidate events, i is the identifier of the +th event, and 
Pk are the yield and probability density function (PDF) of component k, which corresponds 
to the signal and each aforementioned background. 

The signal PDF is determined using MC simulation and control-sample data. A Gaussian 
is used as the Mbc PDF. Since AMb is correlated with Mk-k but not with AMb — Mk-k, 
we use the product of the M^n and AMb — PDFs, each of which is the sum of two 
Gaussians. The parameters of the M^n PDF are extracted from the inclusive data 
sample with cms momentum less than 1.5 GeV/c, using the same selection criteria as 
those for B decays. The mean and width of the Mbc PDF and the main Gaussian of the 
AMb pdf (integrated over Mk-k) are calibrated using & —>■ J/'ijjK*^{K*^ —> iP+vr”) 

control sample. 

The B~^ —»• signal PDF has an additional component where a low momentum 

pion is incorrectly assigned as the pion from B decay. We construct its PDF from the 
product of the double-Gaussian in M^n with a two-dimensional (2D) smoothed histogram 
in Mbc and AMb- The fraction of this component is estimated to be (37.4 ± 0.7)% from 
signal MG and is fixed in the fit. 

For nonresonant background, the Mbc and AMb PDFs are taken to be the same as the 
signal PDFs integrated over Mk-k, while a second-order polynomial is used for the Mk-k 
PDF. For combinatorial background, a threshold function [ig is used for the Mbc PDF. For 
B^ —> J/ipD^, a first-order polynomial is used for the AMb PDF. To take into account the 
kinematic boundary for B^ —>• J/'0D°7r+, we use another threshold function. 


Fthr(a^,Xc;p, c) 


{x — XcYe (x > Xc) 

0 {x < Xc) 


with X = AMb and Xc = M^n — (?^_b+ — — m.,r+), where mB+ and are the nominal 

B~^ and masses. The same M^n PDF as used for signal is used for cmb(DO), while a 
first-order polynomial is used for cmb(fake DO). For peaking backgrounds, the PDFs are 
modeled by 3D smoothed histograms from a large J/ip inclusive MG sample. 

In the fit, the value of Nk and the parameters for the polynomials and threshold functions 
are allowed to float. 

Figure HI shows {Mkw-, AMb) scatter plots for candidates in the Mbc signal region. There 
are 11 candidates in the signal region for B^ —>■ J/ipD^ and one for F+ —>■ J/ipD^ir^. 
Table HI summarizes the results. The signal yields (F) in the signal box, the expected 
total background yields (6), and their statistical errors are obtained from the maximum 
likelihood £t. The efficiencies (e) are determined from the signal MG sample with the 
same event selection used for the data. A three-body phase-space model is employed for 
F+ —>• J/ipD^n^ decay. The AMb-, Mxn and Mbc distributions are plotted in Fig. |21 
Also shown are projections of the maximum likelihood £t result, which agree well with the 
data. No significant signals are found. We determine 90% G.L. upper limits for the signal 
yield (Fgo) and branching fraction (B) from the observed number of candidates (no) and the 
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FIG. 1: {Mktt, ^Mb) scatter plots for data in the Mbc signal region for (left) and 

—> J/ipD^TT'^ (right). Dashed lines indicate the signal regions for Mk-k and AMb- 


estimated backgrounds (6) in the signal box using the Feldman-Cousins method Q. The 
systematic errors due to the uncertainties of signal detection efficiency and background yield, 
elucidated below, are taken into account Q. The decay branching fractions i3(J/'i/i — > ) 

and B{D^ —> K~7i^) are taken from the world averages 3|. The fractions of neutral and 
charged B mesons produced in T(4S') decays are assumed to be equal. 


TABLE I: Summary of results: Y and b are the signal and expected total background yields in the 
signal box, no is the observed number of candidates in the signal box, e is the detection efficiency, 
Ygo and B are the 90% C.L. upper limits for the signal yield and branching fraction. 


Mode 

Y h 

no 

6(%) 

Tqo 

B(10-5) 


-l.OtJ)' 14.6±1.2±0.6 

11 

29.9±2.8 < 4.0 

< 2.0 

5+ 

2.36±0.36±0.21 

1 

14 Q+2-8 

< 2.5 

< 2.5 


The systematic errors on the background yields are evaluated by varying each hxed 
PDF parameter by ila of the measured error, by increasing the order of polynomial for 
combinatorial background, and by changing the Mk-w PDF of nonresonant background to an 
exponential function. The changes in the background yields induced by individual variations 
are added in quadrature. The systematic errors for the efficiency (Table m consist of the 
uncertainties in tracking efficiency of 4.0% for and 6.1% for B^ mode, in particle and 
lepton identihcation of 2% per track, in branching fractions B{D^ of 2.4% and 

B{J/'ip —>• of 1.7%, and in MC statistics of 1.2%. For —>• J/'ipD^n'^ decay, an 

additional systematic error of -|-13.7%/-38.3% due to the three-body phase-space model is 
assigned to the maximum efficiency variation among the slices of M{J/ip, D^), M{J/'ip, 
and M{D^, 7r+). 

We also obtain the branching fractions of the corresponding nonresonant decay channels 
in the 1.71 < Mk+tt- < 2.01 GeV/c^ region from the yields of the nonresonant components 
in the ht. The yields in AMb and Mbc signal region are SO.Olg^s^ for B^ —»• J/'iIjK^'k~ 
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FIG. 2: AMb, Mk-k and M^c distributions in all or related variables’ signal regions, as labeled 
inside each plot, for —> J/'ipD^ (top) and B~^ —> (bottom). Superimposed on the 

data are projections of the signal and summed background components of the maximum likelihood 
fit result. 


and lO.l^lg for —> J/xI)K^'k 7r+. The efficiencies are verified to be the same as for 

the modes with an intermediate resonance. The systematic errors are estimated by 
the same procednre. In decay, we snbtract the contribntion from B^ —>■ X(3872)i^’'' 
[X(3872) —>• J/'^TT+TT”], which is estimated to be 1.20 ± 0.33 candidates by MC simnlation 
with the branching fraction taken from Ref. Q. Finally, we obtain B{B^ —>• J/ipK^n ) = 
(1.51i(o;(g ± 0.15) X 10“® and B{B^ — > J/ipK^T:~TT^) = (O.OlJ;® ± 1.6) x 10“® in the limited 
Mk+tt- region (where the hrst errors are stat. and the second are syst.). 
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TABLE II: Summary of the contributions to the systematic uncertainty (%) on the detection 
efficiency. 


Source 

RO 

B+ 

PID efficiency 

8.0 

10.0 

Tracking efficiency 

4.0 

6.1 

MC statistics 

1.2 

1.2 

J/ip branching fractions 1.7 

1.7 

branching fraction 

2.4 

2.4 

3-body decay model 

- -M3.7/-38.3 

Total 

9.5 

40.2 


In summary, we have performed a search for > J/'ipD^ and —>■ J decays. 
No signal is observed for either decay mode and upper limits on the branching fraction at 
90% C.L. are determined to be 


B{B^ J/^pD^) < 2.0 X 10-^ 

B{B+ < 2.5 X 10“^ 

The results are consistent with the BaBar results and rule out the explanation of the 
excess in the low momentum region of the inclusive J/V' spectrum as intrinsic charm content 
at the 1% level in the B meson. 
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